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^ ■ Abstract 

oo ; 

CN \ It is shown that in the presence of a background magnetic field within solar 

' interior a cavity for low frequency MHD eigen modes (with periods 1-10 days) 

J> ■ near equatorial plane can arise. The lower boundary of the cavity coincides with 

ON , the center of the Sun while the upper one corresponds to the Alfven resonant 

| layer where high accumulation of wave energy takes place. The localization z s 

CN) ■ and the width of the Alfven resonance layer Az s are determined by: (i) the node 

number of eigen modes n = 1,2,..., (ii) by the angle a of oblique propagation 
of waves with respect to the magnetic field B, and (hi) by a low magnitude of 
the background magnetic field itself, Bq = 1 — 100 G. The amplitude of eigen 
O { oscillations in a resonant layer determines the density fluctuation value dp/p 

^ | that is restricted through the imaginary part of eigen frequences. 

i— t ' For large node numbers n>l there appear many narrow resonant layers 

where a neutrino propagates through a large density fluctuation 5p/p with 
the oscillation length that is much bigger than the width of a resonant layer, 
. ^ ' l osc = AirE/Am 2 3> Az s (n). It is shown that neutrino crosses many such 

^ ■ bumps on the exponential background profile that motivates to consider these 

MHD waves as a plausible matter noise for the MSW solution to the Solar 
Neutrino Problem (SNP). 
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1 Introduction 



The existence of the solar neutrino problem is now established from five underground 
experiments 0, [| f|, [J more or less independently of any details of the Standard 
Solar Model (SSM) @. 

At the same time it becomes clear ]7| that any modifying SSM astrophysical so- 
lution contradicts the reconciled data of the Homestake |J and the Kamiokande 
experiments. Moreover, Helioseismology data are in a good agreement || with astro- 
physical parameter values given by SSM and relevant to all Particle Physics solutions: 
(i) neutrino vacuum oscillations ||, (ii) resonant neutrino oscillations in medium or 



the Mikheev-Smirnov-Wolfenstein (MSW) solution [K| and (iii) the Resonant Spin- 
Flavor Precession (RSFP) in medium fIT||. 

The most popular Small Mixing Angle (SMA) MSW solution|T(J depends on min- 
imal number of unknown fundamental parameters (neutrino mass difference Am 2 and 
mixing sin 2 29 <C 1) and relies on resonant neutrino conversions [y^ — > v^l) i n matter 
with the smooth density profile p = p 00 exp(— r/H) where H = 0.1R Q is the height 



scale in SSM |18||. Any change of this profile generally may be crucial for allowed fun- 



damental parameter region (see below) while standard Helioseismology corrections to 



SSM give density fluctuations deep in the solar interior at the level Sp/po < 1 %|12 
that occurs too small to influence neutrino oscillations. 

Another scenario (RSFP JTT|] ) includes the additional parameter pBo where Bo is 
an extrapolation of a regular large-scale magnetic field (as seen on the solar surface) 
down to the bottom of the convective zone and p is an unknown transition Majorana 
(Dirac) neutrino magnetic moment. The existence of such large-scale magnetic fields 
Bq under the bottom of the convective zone and even in the solar core is currently 
discussed (see, for instance, |13|), and it is not excluded that these fields can lead 
to the MHD wave generation in the solar interior relevant to the MSW-resonance 
region for neutrino oscillations since density perturbations Sp(f,t)/p ( as well as the 
magnetic field ones, b(f,t)/B , and the speed fluctuations v(f,t)) are evolved in the 
presence of a central field B . 

Therefore, the main goal of present work is the search of possible magnetic field 
enhancement of the density perturbations 5p/po resulting in an effective non-direct 
neutrino interaction with an external magnetic field B Q without any neutrino magnetic 
moment. R 



3 Another example of neutrino interaction with magnetic field and also without neutrino magnetic 
moment was found in fl4|| where the mean axial vector potential Va ~ GF^elzlb^) ~ Bq for an 
electron gas polarized in a strong magnetic field B = (0, 0, Bq) jL5j was calculated. This polarization 
leads to anisotropic neutrino emission through MSW oscillations studied in fl6|| in order to explain 
kick of a neutron star. For the solar case the axial vector potential is not efficient since the solar 
magnetic field is too small to polarize the Boltzman electron gas. Even at the bottom of the 
convective zone with the strongest magnetic field value Bq ~ 10 5 G the paramagnetic contribution 
is [LbBq < 10~ 3 eV only and the axial vector potential Va is much less than the vector potential V, 
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Here we consider only very modest values of large-scale background magnetic fields 
in central region of the Sun, Bq ~ 1 — 10 — 100 G, that is not in contradiction with 
a primeval origin of such central field in the stars of main sequence like the Sun and 
with non-observation of traces of such small (dipole) fields in the solar photosphere. 
This is in contrast with the visible traces of more stronger toroidal (and poloidal) 
fields generated at the bottom of the convective zone. 

In particular, the most interesting region B < 10 G where MHD density pertur- 
bations are similar to the intense matter noise (see below) obeys the constraint on the 
central field Bq < 30 G found recently |17| from the assumption of primeval origin of 



solar magnetic fields which are evolved via dynamo mechanism. 

There were some previous attempts to modify the SSM density profile po JH 
entering the neutrino evolution equation in the MSW solution to SNP. The influence 
of periodic matter density perturbations above the average density po, i.e. 

p(r) = p [l + hsm(n/r)] , (1) 



on resonant neutrino conversion was investigated in Ref. f[5 |. In this case the para- 
metric resonance effect appears when the fixed frequency (7) of the perturbation is 
close to the neutrino oscillation frequency and a large amplitude value (h ~ 0.1 —0.2) 
is assumed. There are also a number of works which address to similar effects by 
different approaches PD[ . 

One of disadvantages of models like Eq. (|]) is our unawareness of values of model 
parameters 7, h as well as of physical origin of these perturbations. 

Another approach was developed in Ref.|^T[ where authors assumed the presence 
of a white noise of matter density ((£) = (5p/po) = 0) added to po instead of the 



second term in Eq. (0). The number of free parameters (the noise amplitude J (£ 2 ) 
and the correlation length Lq for random 5-correlated domains) remains the same one 
(=2) while the presence of white noise in the Sun is doubtless since there are many 
mechanisms to enhance random density perturbations. 

However, Bamert et al (c.f.,for instance, [p2| |) stated that in analysis of matter 



noise influence the MSW solution to SNP given in |2T| (as well as in their own first 



work and in [23]) what has been missing so far is a plausible source for such 5- 
correlated fluctuations in the vicinity of the MSW resonant point. 

Assuming that a certain superposition of g-modes can be a source of such fluctu- 
ations deep in the solar interior, they showed that even for large amplitudes of some 
g-modes (or large Sp/p ~ 4%) for a large size of density fluctuation cell, L > 3000 km, 
the direct numerical integration of evolution equation leads to standard MSW result. 

Comparing this value with the neutrino oscillation length at the MSW resonance, 
KVc = lv/ sin 2$ = 250 km(E / MeV) / Am\/ sin 26, they found that the main assump- 



tion of |21] about the (^-correlated form of matter density fluctuations, L <C l 



Va ~ VplbBq/T <C V where the standard MSW vector potential is defined as V — GFV2pY e /m p 
for v e — > v p i.e. it is proportional to the mean matter density p. Here Gf is the Fermi constant, \xb 
is the Bohr magneton, T is the temperature within the Sun, Y e is the electron abundance. 



3 



is violated. Really, the opposite inequality holds, L ^> l osc , especially for typical 
g-mode wave lengths, X g ~ L > 0.1R Q = 7 x 10 4 km if the scale of random domains 
is associated with a g-mode wave length. 

Let us turn to the MHD waves trapping in cavities which are appearing in the 
central regions of the Sun as a result of simultaneous acting of magnetic fields and 
density inhomogeneities. There are two advantages of these oscillations: (i) accu- 
mulation of wave energy in a resonance Alfven layer that leads to a large Sp/po (see 
below); (ii) short wave lengths Xmhd <^ lose that allows us to consider these waves 
as a source of 5-correlated matter noise. 

First, notice that strictly speaking these MHD waves are called magneto- acoustic- 
gravity waves (MAG) |24| (see old bibliography therein). The spectrum of MAG 
waves is very wide and near the photosphere some of them can be trapped in the sub- 
photospheric waveguide with consequent transformation into the magnetosonic waves 
characterized by 3-5-min periods which indeed are observed in the active regions of 
the solar atmosphere. 

We are interesting here in low-frequency branches of MAG waves which can lead 
to long periods (like days or weeks) for temporal variations of neutrino fluxes that, 
in general, can be observed in the SuperKamiokande (SK) experiment. 

In the recent paper by Guzzo et al||25[l it was found that in the presence of gravity 
force ideal magnetohydrodynamics (MHD) has a solution for magnetic field pertur- 
bation with the periods of a few days that correspond to slow and Alfven waves, 
ojs ~ — kB Q / ^Anp{r) where k is the wave number (free parameter obeying 
k > 2tt/R q , Rq is the solar radius), B is a regular magnetic field in the Sun and 
p(r) is the SSM matter density. 

A search for such periodicity in SK could be, in principle, realized. However, the 
authors of[^| confined themselves with the statement that "...time fluctuations of the 
solar neutrino flux have to be observed if the resonant neutrino spin-flavor precession 
solution to the solar neutrino problem works" . 

In contrast to PS[, we show that for the same MAG waves while one neglects 



neutrino magnetic moment or RSFP [[11]] is absent the MSW solution to SNP is 
nevertheless crucially modified and can lead to the corresponding temporal variations 
of neutrino flux. 

In addition, the authors [^ did not calculate in a self-consistent way the location 
r s and the width Ar s of the singular (resonant) layer since they did not solve the 
eigen value problem on the base of solutions of the Heun equation (Eq. (|15D below) 
or they did not find the corresponding MAG wave spectrum (including the imaginary 
part of frequency connected with Ar s ). 

We show that low frequences of eigen oscillations correspond to the resonant layer 
position deep in the solar interior (under the bottom of the convective zone) or a 
cavity for MAG waves being bounded by: (i) the center of the Sun and (ii) the Alfven 
resonant layer is localized in the central region where matter density oscillations can 
influence the MSW resonance for neutrinos. 
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Note that large density fluctuations Sp/po in these narrow layers, Az s <C l osc , can 



influence the MSW conversions like a matter white noise [21] if these layers are tightly 
placed along neutrino trajectory. 
This paper is organized as follows. 

In section 2 we present the total set of MAG equations and derive the master 
Heun equation. 

In section 3 we find asymptotic solutions of master equation near the center (sub- 
section 3.1), at the surface of the Sun (subsection 3.2) and around the singular Alfven 
layer (subsection 3.3). 

In subsection 3.4.1 we discuss the cavity model in ideal MHD and collisionless 
damping of MAG waves in an Alfven resonant layer. Then in subsections 3.4.2, 3.4.3 
using reasonable boundary conditions from the common solution in whole region 
< r < R Q we derive the dispersion relation for the discrete frequency spectrum 
cu n (B , k x , k y ). Here the order of modes n = 1,2,.. enumerates the decomposition 
over the set of MAG eigen waves replacing Eq. ([!]) to 

p{z,t) = p {z)[l + Re J28 Pn (z,t)/p (z)] . (2) 

n 

Then in section 4 we calculate the resonance layer location z s and its width Az s (n) 
making use of the wave energy flux bounded through imaginary part of eigen fre- 
quency. 

In section 5 we present main results for matter density fluctuations that should 
be substituted into Eq. (||) and calculate the enhancement of a density perturbation 
Sp n /pQ in a resonant Alfven layer. 

In final section 6 we discuss the validity of our results and state some problems 
for future exploration. 



2 Master equation for MAG waves 

For MAG waves we use the full set of MHD equations including the gravity force 
pg where g = V$ and the gravity potential $ is given by the Poisson equation 
A0 = —AirGp, G ~ Mpf is the Newton constant. 

The Euler equation (conservation of momentum) is of the form 

dv 1 r -> 

P dt =-VP + P9 + ^[ rotBxB \ > ( 3 ) 

where d/dt = d/dt+v-'V is the substantial derivative and the last (Lorentz force) term 
may be transformed to the sum of the magnetic field pressure term —'VB 2 /8n and a 
perturbation magnetic field contribution — (47r) -1 f[& x [V x B]] + [Bx[Vx b]]j when 

in linear MHD we decompose the total field B = Bq + b(t, r) (see below subsection 
2.1). 
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The mass conservation leads to the continuity equation 

dp 



dt 



+ pu = 



(4) 



Here for compressible gas the function u = div v 7^ is the main object for search of 
solutions with MAG waves. We decompose density in linear MHD as p = po + Sp(t, r) 
with the background density p given by SSM [|IBf . 
The energy conservation obeys 



dp pdp 

Tt- 1 - p Tt=- {l ~ 1)Q 



(5) 



where the factor 7 = c p J Cy is given by the heat capacities, Q is the rate of all energy 
density sources and losses. 

Finally system is completed by the Faradey equation for the magnetic field, 

dB -> c 2 -» 

_ = rot \v x B] H A£? , 

9t L J 4tto- 

that transits in ideal MHD (the conductivity a — > 00) to the equation for frozen-in 
magnetic field, 



db „ 



v x B 



(6) 



In the last equation we assumed that the background field _B does not depend on 
time. 

Obviously, these equations transit to the Standard Helioseismology Model (SHM 
) equations in the limit B —>■ 0. 

The system of linearized Eqs. for non-adiabatic perturbations in a com- 

pressible ideal conducting medium with arbitrary direction of the magnetic field B 
and in the Cowling approximation (neglecting the perturbation of the gravitational 
potential, Scj) = 0) is reduced to 



d5p 



+ (vV)po + p u = 



dt 

dv 1 r -> -* -» -• 

po di + V6p ~^ 6p+ A^ i v(jB ° 6) " (jBoV)6 

( < ' +(vV)T ] + Po u = 5Q , 



dt 



PoCv 

^ - (B V)v + B u = 
5p _ 5p 5T 
Po Po T 



(7) 
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Here v, 5T, Sp, Sp and b are small Euler perturbations of the velocity, the temperature, 
the density, the pressure and the magnetic field strength correspondingly. One sup- 
poses that at initial state the gas is immovable (vq = 0) and for the constant magnetic 
field B the ideal gas (p = Rp T Q ) obeys the hydrostatic equilibrium, Vpo = p$g. 
Notice that fifth equation in follows from equation of state of ideal gas. 

In what follows we consider the wave dynamics in central regions of the Sun 
where the adiabatic condition, SQ = 0, is fulfilled with good accuracy for the high 
frequency oscillations in SHM and MHD. To simplify the further analysis we adopt 
this condition for whole spectrum of oscillations in MHD. 

Since we assume the background magnetic field B has the dipole structure one 
expects that it occurs both constant and uniform near the dipole axis or we can put 
B = const within central region of the Sun where a cavity for eigen MAG oscillations 
is localized. In addition to, since neutrinos propagate to the Earth along the ecliptic 
one can suppose that the dipole magnetic field is perpendicular to the direction of 
neutrino momentum. 

Thus, we choose the Cartesian coordinate system with z-axis directed upon op- 
posite direction to the neutrino momentum or along g. This means: (i) z — at the 
surface of the Sun and (ii) z = R Q at its center. In this system the magnetic field 
and acceleration vectors are B = (B o ,0, 0), g = (0,0, g(z)). 

In order to obtain an analytical solution of the system Eq. (0) we adopt addi- 
tional assumptions based on an effective short wave length along z-axis, X z <C H p —\ 
(1/ po)dpo/dz | _1 where H p = H = 0.1R Q = const is the density height scale fll8 



Really in most interesting cases of cavities in deep solar interior the MHD wave spec- 
trum u n depends on large order (node) numbers n ^> 10 with the effective wave 
length ~ L/n where L ~ z s < R Q is the cavity width, z s is the Alfven resonant 
layer position (see below section 4). 

On the other hand, using the real model of the Sun [^6| one can easily show 
that both the temperature T (z) and the acceleration g(z) have large height scales, 
H T =| (l/T )dT (z)/dz I" 1 and H g =\ (1/ ' g)dg(z) j ' dz | _1 correspondingly, obeying the 
inequalities H g , Ht > H p . As result we may consider the gravity force acceleration g 
and the temperature To as some constants on scales of MHD wave lengths, g ~ const, 
Tq m const. This approximation remains valid even on whole cavity width L for 
"narrow" cavities 3H < L < H T , H g including the MSW resonant point r MSW ~ 
O.3-R in the solar interior while it fails near the center of the Sun, z -C H p , where 
9 ~ z. 

The isothermal approximation with g pa const means that the density and the 
pressure are given by the exponential (Boltzman) law: po ~ e z / H , po ~ c z l H obeying 
the ideal gas equation state po/po = RT Q / p with the gas constant R = K,/m u and 
the molecular weight p measured in mass units m u ~ m p . Here k is the Boltzman 

constant; H = c 2 /g = const is the height scale fllBfl; c = yRTo/jp = const is the 
isothermal sound velocity and the adiabatic sound velocity is given by c\ = 7c 2 = 
const. 
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Thus we are looking for solutions in the form v = v(z) exp(^—i(ujt — k x x — k y y)^j 

(and analogously for 5p and b). We keep here notations of the paper f2T| . In these 
notations we obtain from the system (|7J) all perturbation functions. In particular, 
velocity perturbations are of the form 



ik x (c 2 s u + gv z ) 



uj 



2 , ,2 



v A uo 



uj 2 — k 2 v\ 



u 



(9) 



( 



uj 2 g 



k\g 2 - cj 4 



i _ k 2 ± c 2 , - 1) - ky A 



x{ 



+ 



uj 2 v\ 



uj 2 — k 2 v\ 
1 du~ 



u + 



klvV gdz 



■}• 



(10) 



^From the continuity equation Eq. (Q) or from the first equation in (0) we have 
obtained the density perturbation, 



5_p 
A) 



1 



i ^ 



and from Maxwell equations for ideal plasma Eq. (^) or from the fourth equation in 
(|7J) we find magnetic field perturbations, 



b x 



-DO + —-DO ? 

UJ too 

k x v y 
'--do , 



UJ 



B n . 



(12) 

(13) 
(14) 



One can easily see from Eqs. (p|)-(|T4|) that all perturbation functions are expressed 
via one unknown function u(z) = dv z /dz + ik x v x + ik y v y . 

The master equation for u(z) follows from last definition and coincides with Eq. 



(3) in [27 1 



+{£[(£ - l)b + K{[a + —-) - lj - K{a}u = 

I X 



(15) 



More than hundred years this equation for an arbitrary argument £ is known in mathematics 



as the Heun equation |28 
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where £ = ur jk x v A is the independent variable; Va(z) = B / JAttpq(z) is the Alfven 



velocity; K± = yK 2 + K 2 and u are the transversal wave number and the oscillation 
frequency correspondingly. 

The dimensionless parameters in Eq. (|T5|) are given by 



1 0= + K, 



yK 2 > 7 " V 7^ 2 

ii — — — , K x — k x H, Ky — kyH, 

£ = -W^)5 fa) = 4 = Poe z/H . (16) 

7 A x W A 

Eq. (|T5|) is the ordinary differential equation of the second order with four singular 
points: £ = 0, 1, a, oo. The solutions of this equations are given by the Heun 
functions P8|. These solutions were explored in the paper however, they are too 
complicated for application to solar neutrino physics. 

Therefore, we have simplified Eq. fll5|) in order to find solutions through elemen- 
tary functions and to study main physical properties of the MAG wave cavity within 
solar interior. 

In general, when both K x ^ and K y ^ 0, the master equation Eq. (|T5D has two 
singular levels: the resonance level £ = a arises for the condition oj 2 = c 2 s v 2 A k x / {c 2 +v\) 
(casp resonance) and the condition £ = 1 is fulfilled for uo 2 = k 2 v\ {Alfven resonance) . 

As shown in |27j for the solar conditions the casp resonance level lies in medium 
above the Alfven level, however, the distance between these levels, AzAcasp — Hln[l — 
Q 2 /^fK 2 ]~ l is negligible for low frequencies Q -C because we are interesting in 
1-10 days variations of density fluctuations with the probable influence neutrino flux. 

When waves reach resonance level which has a finite width in z-direction and 
is infinite in both perpendicular directions they can be efficiently absorbed within 
resonant layer. This leads to concentration of wave energy along magnetic field. The 



energy absorption coefficient is of the form [24] 



D = 1 - exp(-2 7 r v / 46- 1) , (17) 

or absorption becomes more efficient for low frequencies, Q — > 0, b — ► oo, D — > 1. 

The coincidence of casp resonant layers with the Alfvenic ones for low frequences 
(il lK 2 ) means that the resonant layer £ = 1 may appear only within the Sun 
where the condition va c s is fulfilled that corresponds to the inequality for gas and 
magnetic field pressures, p gas = p ^> p ma g = B 2 /8n. 

In such conditions we can simplify our master equation Eq. ([15]) that can be 
written in the final form 



5 In this formula one neglects the temporal decay of a wave amplitude since free oscillations 
without lower reflector and thereby with zero imaginary part of the frequency Im f2 = are assumed. 



9 



+{£[(£- l)b + Kl{\ + 
Here one supposes also 7 ~ const. 



- 1 



#|}« = . 



3 Asymptotic solutions 

The main idea how to find simple asymptotic solutions of Eq. ( jOf) comes from its 
physical simplifications in three important regions: near the photosphere z = 0, 
I £ |<C 1, around the singular (Alfven) layer | £ — 1 |< 1 and near the center of the 
Sun, z = Rq where | £ |^> 1 for the same fixed frequency u. In these regions solutions 
are expressed through elementary functions. Using conditions of linear dependence of 
these regional solutions we build a more general solution that is valid in whole region 
—00 < z < 00. 

Then using boundary conditions at the center of the Sun and at the solar surface 
we have found eigen frequencies (subsection 3.4). 



3.1 Solution near center of the Sun 



In central region of the Sun for | £ |> 1 including the center, z s < z < Rq, where the 
resonant Alfven layer position z s is given by the value (u)/k x VA{z s )) 2 = £ = 1 we can 
reduce Eq. (|H|) to the more simpler form 



General solution of this equation is of the form 

1 

71 



u = Ul = ^ 7 =(a 1 C u + A 2 C tu ) , (20) 



where v = \jb — 1/4 w (K±/£l)\J (7 — lj/7, A\p are arbitrary constants which are 
determined from boundary conditions. For low frequencies | v 1 the derivative is 
substituted by 

Since the variable £ is proportional to the SSM density profile £ ~ e z ^ H [p8l the 



solution Eq. (POf) describes the incident wave and the reflected one with the wave 
number k z = v/H. 

One can easily see that the asymptotic Eq. (|I9| ) is similar with the SHM equation 
for g-modes where the wave number b/H 2 w k 2 g("f — Vj/^Hoj 2 in Eq. (ITBp includes 



the buoyancy frequency term N 2 / u 2 ~ ^(7 — 1)/ Huo 2 coming in SHM equations. This 
corresponds to the limit Bq — > when the variable £ increases (£ — > 00). 
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3.2 Solution in vicinity of a resonant Alfven layer 



For £ — > 1 Eq. (WW) transits to 

+ 



(£-1)11^)^ + ^1-1 



u = 



(22) 



Two linear independent solutions of this equation are expressed through the Bessel 
functions of the 1-st kind: 



u 



u 2 = C 1 J fM {^) + C 2 J^) 



(23) 



2u VT 



where we denoted the argument and the index of Bessel functions as \I/ 
jx = 2[1 — If^yl^] 1 / 2 correspondingly and (7^2 are arbitrary constants. 

For the case of MAG wave propagation along magnetic field (k y = 0) the index 
fj, is real, li 2 = 4(7 — l)/7, while for the particular case of oblique wave propagation 
k y /k x = "1/I57/I6 — 1 = 3/4 (ideal H-atomic gas, 7 = 5/3) the real index li = 1/2 



(24) 



^157/16 — 1 = 3/4 (ideal H-atomic gas, 7 = 5/3) the real index li 
corresponds to the trigonometrical functions instead of the Bessel ones, 




Ci sin * + C 2 cos * 



As it follows from the decomposition of Bessel functions near the resonance level 
£ ps 1, or explicitly for the particular case ( |24] ) the second independent solution 
diverges in the limit £ — ► 1. However, since this singular behaviour is an intrinsic 
property of the solution Eq. Q2HD caused by physical reasons we may not exclude this 
singularity by the forced condition C 2 = 0. It is important for our problem to have 
convergence (regular solution) in all points along integration path | z |— > 00. An 
analytic continuation (transition to the complex argument £ = £1 + ^£2) originated by 
a physical mechanism discussed below means that £ = 1 is the removable singularity. 

Notice that the solution Eq.([23|) is oscillating along the direction to the center, 
£ > 1, and has exponential behaviour in opposite direction (to the surface of the 
Sun), £ < 1. First derivative of the solution Eq. ( p3[ ) is of the form 



du 2 2h> 2 
where fJH?) = (li^)J^) - J^V). 



C.j'm + CjA^) 



(25) 
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3.3 Solution near surface of the Sun 



The argument £ = uj 2 / (k 2 v\) in Eq. ( [L6|) is given by the parameter f3 = Poe z ^ H where 
near the solar surface z = the factor (3q is changing in a wide region] 



40, 
{ 400, 
40000, 



for 
for 
for 



B = 
B -- 
B 



1000 G, 
■ 100 G, 
= 10 G. 



(26) 



However, for any /3q and K x > 1 the argument £ occurs small near photosphere, 
| £ |<^C 1, since we have considered low frequencies through whole solar interior, 
uj < 10 -5 As result the master equation Eq. ( |18|) transits near photosphere to 



^ 9 d 2 u an i 



du 



d£ 



. 7 -IK* 



n 2 



7 



+ K\ 



u 







(27) 



which has the solution 



u 



U 3 = DiJ 2 K ± (y) + D 2 J-2K ± (y) 



(2i 



Here the Bessel functions depend on the argument y = Iv^J— £, D\ 2 are new integra- 



tion constants. The demand of regularity of the solution Eq. (28) on the photosphere, 
£ — ► 0, results in D 2 = 0. The solution Eq. (E^) describes evanescent oscillations 
with an exponential decreasing amplitude at the surface of the Sun. 
The derivative of the solution Eq. ( p8|) is given by 



duz 



-Di 



2u 2 



2K\ 



J2kAv) - J 1+2kAv) 



(29) 



3.4 Spectrum of MAG waves 
3.4.1 Cavity model in ideal MHD 

We have shown that the Alfven resonant level £ = 1 divides whole wave propagation 
space into two parts in which MHD waves have distinct properties. In the region 
1 < £ < +oo the solutions u%(^) and w 2 (£) describe waves that propagate in both 
directions along z-axis while in the region < £ < 1 the solutions U2(0 and 
describe evanescent oscillations that decrease in an exponential way to the surface of 
the Sun, z — 0. The level £ = 1 itself is the singular level where %(£) — * oo. 

Therefore, waves which are propagating from the center of the Sun reach the 
resonance level £ = 1 where their partition implies reflection back and partial capture 
at the resonance layer. In general, high wave amplitudes at the resonance imply a 
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start-up of different dissipation mechanisms (linear or non-linear) within that layer. 
The condition v z = for the hard center of the Sun is a natural demand due to high 
pressure there and leads to full wave reflection from the center resulting in a cavity 
appearance with the lower boundary at the center z = R Q and the upper boundary 
z = z s at the Alfven resonance layer. 

The cavity is defined within the region 1 < £ < ^ and its width L = Hln | ^ | 
is given by the argument value at the center of the Sun, ^ = (£l 2 /'jK^)P e R@ ^ H . 
Notice that the width L depends significantly on the value of a central magnetic field 
via the parameter /?o in Eq. (^). 

The waves which are trapped within cavity loose wave energy due to partial ab- 
sorption of their energy in the resonance Alfven layer with the rate given by an 
imaginary part of the eigen frequency Im (u) ^ 0. 

This dissipationless (i.e. reversible) damping is analogous to the collisionless mech- 
anism of coherence losses for eigen (Alfvenic) modes considered in |30| for the Alfven 
resonant layer which is the sheath of a solar coronal loop. The resonant absorption of 
Alfvenic surface waves leads to the heating of solar coronal loops. Moreover, the iner- 



tia of electrons m e ^ oo (v e x = yT/m e ^ 0) should be taken into account or kinetic 
approach for wave-particle interaction is necessary to explain collisionless damping. 

The same collisionless damping of MHD waves with absorption of waves in the 
Alfven singular layer as a mechanism of heating of laboratory plasma was considered 
in |3l| and in [[3^] - for heating of a nonuniform cosmic plasma. 

In ideal MHD for a non-uniform medium one can also obtain the collisionless 
damping of eigen modes from the dispersion relation solving the corresponding eigen 
value problem with given boundary conditions [30]. 

There is another analogy that is well-known from plasma physics where the col- 
lisionless decrement of plasma waves comes both from dispersion relation for eigen 
modes and from the kinetic theory beyond the ideal MHD. The collisionless damp- 
ing and increment are equal to each other for rates of the direct (inverse) Cerenkov 
emission (absorption) e «-> e7 in the wave-particle interaction. 

The Landau damping of the Alfven waves u = k x VA(z s ) in the inverse Cerenkov 
process , 7^e — > e, starts only within the Alfven resonant layer £ = 1 obeying the 
standard (Jerenkov condition uo = k ■ v e where the electron velocity v e is much bigger 
than the Alfven one, v e ~ v e T 3> v a- This means that those electrons which are 
moving almost across Alfven layer (in ecliptic) could absorb Alfven waves . 

Thus, accounting for trapping of waves in the cavity 1 < £ < £oo with partial losses 
of the wave energy within its upper boundary (£ = 1) one can solve the problem of 
infinite wave amplitude at the resonance singular layer in framework of ideal MHD 
without taking into account of usual dissipation mechanisms (viscosity through colli- 
sions, ohmic dissipation, etc.). Really, for complex frequency the argument £ becomes 
complex or the singularity £ = 1 is removable (see below). 

Our main task in present work is the derivation of complex eigen frequencies and 
eigen functions in a cavity we have considered. First step for that is the construction 
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of a common solution in the whole variable region | £ |< oo. 



3.4.2 Common solution and continuity conditions 



Let us assume that the complex frequency Q — Q± + iQ 2 = + id) obeys the 
stability condition Vt 2 < § (Q\ > Q 2 ) and input the following notations for the 
independent variable £ = £1 + ^£2 = £i(l + ^C) : 



6 



(5{z) 



6 



^1 ^ 



nj 

2d 



(l-d 2 )P(z) 



l-d? 



(30) 



Notice that transition to free oscillations corresponds to the limit d = Q 2 /Qi 0. 

In the Fig. 1 we plot definition region for all three solutions upon the axis £1 = Re £ 
where the region < £1 < £100 is the whole integration region. In common regions 
for overlap of solutions 1*1,2,3(0 we can use analytic continuation and connect these 
solutions each other. 

By continuity of functions 1*1,2,3 and of their derivatives dui^^/d^ we choose two 
matching points in two overlap regions : (i) £1 = 1/2 for the pair of solutions 1*3 (£) 
and 1*2 (0 and (ii) £1 = 3/2 for the second pair u±(0 and 1*2 (0- Then at these points 
we obtain the continuity conditions 



MO = MO; ^ = ^ for £1 



d£ 



fc \ fc \ d M0 du 2 (0 f c 
For the complex matching points chosen above, 

£o = ^(l + <), £* = -(! + <) , 



1 

2 ' 
3 

2 ' 



(31) 



(32) 



the arguments of functions in Eqs. 



are of the form 



*o = *|&=i/2 = 2z/y£ - 1, 2/o = 2/l6=i/2 



6=3/2 



2z/J36 



1 . 



(33) 



Substituting the solutions 1*1,2,3 and their derivatives u 12S given by Eqs. (P0|), (]2T[) , ( [23 
fl25|), (|28|) ,(|29|) into the continuity conditions Eq. (|31|) we obtain finally the system 
of algebraic transcendental equations 



CiJ^o) + C 2 J_ M (* ) = DiJzkAu 



±\yo) , 
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cvj;(*o) +c 2 /_ M (*o) = -^i— 4rx(ito) , 

d J M (^) + c 2 J_ M (*.) = ttp (ACf + a 2 c 

CiJ;(**) + C 2 Jl M (*0 = -^(AiC - AaC*) , (34) 

from which we can define all arbitrary constants except of one, for instance, we choose 
D\ 7^ 0. Let us recall that in the solution near surface of the Sun in Eq. ( Egp we used 
the boundary condition D 2 = excluding irregular part of solution. 

Notice that in low frequency approximation Q <C 1 the index v is large, | v 1 
(see Eq. (pOD), resulting in large arguments of Bessel functions: | |>| \l/o |3> 1, 
| y |> 1. As result we substitute Bessel functions by their asymptotical expansions 
for large arguments. 

3.4.3 Dispersion equation and spectrum of MAG waves 

Let us turn to the boundary condition at the center of the Sun, v z \ z= r q = 0. From 
the Eq. (|T0|) in the limit | £ |^> 1 we obtain the ratio for coefficients Ai )2 in the 
solution Eq. ( POD for central region of the Sun: 

^ = -a- . (35) 

Thus, using the last ratio and the system Eq. (0), we obtain the dispersion equation 
for the MAG waves in the Sun: 



tt* (£-1) 
(e + 1) 



1 , (36) 



where the complex parameter e = S\ + ze 2 = — (£*/£oo) 2 " y can be expressed for small 
damping decrement (d 2 <C 1) via real functions 

e x = -e d * cos $ , e 2 = -e^ sin $ . (37) 

Here the parameter $ = $(fii, k±, a, Bq) depends on the real part of the eigen fre- 
quency Qi. 

3 7 iC 



where v* = Re (v) = \J(l — l)/7(-fG_/^i) is the real part of the index v\ tana = 
K y /K x is different from zero for oblique wave propagation, K y ^ 0. 

As result of the parametrization £1,2(^,^1) in Eq. fl3"7|) the dispersion equation 
Eq. (|36|) can be rewritten as the system of transcendental equations for the real (fii) 
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and the imaginary (d) parts of eigen frequency 

V2 



d(ei-l)+e 2 (l--p) =0 



tte 2 + l-e 1 + -=(l + e 1 ) = 0. (39) 
v2 

One can easily see that d = is not a solution of Eq. (|3£|). This means that full 
reflection from resonance layer never happens. Accounting for notations in Eq. 
we obtain from Eq. fl3"9"| ) another form of dispersion equations 

$ tan $ = — — In 



7 11-6^2' 

d = - In — ^= , (40) 

$ 11-6V2 

where $ is given by Eq. (j38|). 

Since a solution of the first equation in Eq. ( fiO|) depends on the harmonic number 
n = 1,2,... only, $ = $(n) , the damping decrement d depends on that harmonic 
number too, d = d(n). This dependence is plotted in Fig. 2. The calculation gives 
d < for positive values n > 1. One can see from Fig. 2 that number of harmonics 
is finite: 1 < n < n m aa;- While | d | is small for a large n the nonzero value d ^ for 
n max means that energy losses in cavity always happen due to resonance absorption. 

The phase velocity of waves along the magnetic field B normalized to the sound 
velocity V p h/c s = Ui/k x c s = Qi/ '^fiyKx is the second solution of dispersion equations 
Eq. (|40|) as it follows from the definition of the function $ in Eq. (p8[) . This solution 
depends on all set of parameters: the harmonic number n, the value of magnetic field 
i?o, the wave angle of incidence a, tana = k y /k x , the adiabatic parameter 7 and 
the wave number k x . In Fig. 3a we show the dependence V p h(n) on n for different 
values of the magnetic field B Q and for a particular case of the longitudinal wave 
propagation, k \\ B , k y = 0. 

For a given n the dispersion equation has two solutions: the high frequency one 
(plotted by the solid line in Fig. 3a) and the low frequency one (dashed line). The 
cross point of these curves determines the maximal value n max . The value n max 
decreases when the magnetic field B increases. The analogous dependence u>i(n) 
for the fixed magnetic field value 10 G and for different wave propagation angles 
(k±/k x = 1,2,5 corresponds to a = 0°, 60°, 78°) are shown in Fig. 3b. The node 
number n increases with increase of a. 

The change of a node number n in dependence 011 Bo, a is connected with the 
change of a cavity width for these varying parameters. Combining results in Figs. 2,3 
we can find the imaginary part of the eigen frequency f2 = fix + ^2, ^2^) = dQ\, 
for a given harmonic n. 

The peculiarities of the eigen modes ui(n, B , k x , a) are seen from Figs. 3a, 3b: 
the high frequency branch (solid lines) has the minimum value uj^I^ 1 \n max ) for n max 
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and the maximum one ^^al ( n = 1)) for n = 1. Vice versa, the low frequency 
branch (dashed lines) has the minimal value iof^ (n = 1) for n = 1 and the maximal 
one ujj^}(n max ) for n = n max . Note that the first harmonic (n = 1) for the high 
frequency branch corresponds to the magnetosonic waves: V p h ~ c s that is outside of 
our approximation of low frequences, uj -C k x c s . As result we got too fast damping of 
such waves, | d \ > 1CT 2 . 

We do not consider here the low frequency branches (dashed lines) for which the 
rotation of the Sun becomes important. This rotation changes crucially properties 
and spectra of low frequency oscillations (T > T Q ) in the Sun as it was found in P3| . 

Note that since a node number 1 < n < n max is finite a continuum does not arise, 
spectrum is always discrete (ridge structure). 

The spectrum / = LUx(n, B , k x , a = 0)/2tt (in the dimension units fiHz) is shown 
in Fig. 4 for the longitudinal wave propagation: we plotted the real part of frequencies 
u>i(n, B = 10 G, k x ) in dependence on the wave number k x R & for the particular case 
B = 10 G, for k y = and for different harmonics 1 < n < n max = 2780. The 
curve n = 2780 is common for low and high frequency branches. Since the long 
period J" 1 = (1 fiHz)' 1 = 11, 6 days is at the boundary of our approach we conclude 
that all dashed lines in Fig. 4 with n < n max = 2780 and some part of solid lines 
{n > 2000, or / < 1 fiHz) correspond to too low frequencies when we need to take 
into account the Sun rotation (with the period T Q ~ 27 days). 

4 The position z s and the width Az s of a resonant 
layer 

The appearance of a resonant Alfven layer within the Sun, < z s < R Q can lead to 
an interesting application for SNP. From the condition £i = 1 substituted into Eq. 
( pUD we have found the distance of a singular layer from the center of the Sun for the 
case d < 1: 

where Z\ s is the distance of the same singular point measured from the surface at 
z = 0. In Fig. 5 (a,b) we show the harmonic number dependence z s (n) for given 
values of the magnetic field B and the wave incidence angle a. 

One can see that for any values of the parameters B , a and for the low frequency 
branch (dashed curves) the resonant layer is near the center of the Sun (z s / R Q < 0.25). 
For the high frequency branch (solid curve) it occurs at the edge of the solar core 
relevant to the MSW resonance region and either in the radiative zone or in the 
convective one (0.25 < z s /R Q < 1). With increase of the magnetic field B the 
acceptable node number n < n max decreases (see Fig. 5a for fixed k y = a = 0) while 
for the fixed magnetic field value Bq = 10 G it increases if a increases (Fig. 5b). 
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Let us turn to the derivation of the resonant layer width Az s where a significant 
part of the wave energy is concentrated. For this goal we have calculated the density 
of the hydrodynamical energy flux across an Alfven resonant layer 



where the pressure perturbation Sp = (po/iu)(v z /H + ^u) is obtained from the system 
Eq. (0), v* is the complex conjugate value of v z and the function u(z) = 1*2(2) is given 
by the Eq. (£30. 

Notice that we have neglected electromagnetic energy flux density of waves across 
resonant layer because near that Alfven layer (£ ~ 1 ) the z-component of the Poynting 
vector is much smaller than the mechanical one, since the condition 

Pgas 3> Pmag is fulfilled within the Sun. 

In the case of free oscillations the imaginary part of the eigen frequency is zero, 
Cl 2 = d = 0, and the flux pj mec/l )(£) \ s infinite at the singular layer, p( mech )^ = 
1) = 00. However, for eigen modes bounded within cavity (d 7^ 0) this flux is finite, 
p(mech)^ _ -q _^ ^ anc j ^ e function pj mec/l )(£) h as a finite maximum for the 

argument £1 = 1. 

The height of this peak is determined by the value of d = f^/^i- m order to find 
the half-width of that peak we have used the equation P^ mech ' (^i(z)) = P^ mech '(^t = 
l)/2 which has two solutions: £f < 1 and > 1. 

For the half-width we have found 



In Fig. 6 we show the dependence of the width Az s on node numbers n for 
different values of the background magnetic field B and for the longitudinal wave 
propagation, k y = 0. One can see that the width Az s becomes more narrower for large 
node numbers lowering to Az s < 1 km for the case Bo = 10 G in the region n > 1000 
(radiative zone, see Fig. 5a where the corresponding position z s (n) is shown). 

We have found too narrow spikes (hundreds meters) near the node number n ~ 
2680 in the most interesting region z = tmsw ~ O.3i? for the same ambient magnetic 
field B = 10 G while for the field B = 100 G the width reaches Az s ~ 5 km for the 
same position z s with the node number n ~ 268 (compare Figs. 7a, 7b below). 

In reality the position z s (n) is discrete in dependence on the node number n 
(it looks continuous in large-scale Figs. 5a, 5b). For the high frequency branch 
for which our approximation of the Sun at rest seems to be valid we have shown 
in Figs. 7a, 7b (for the background magnetic field B = 10 G and B = 100 G 
correspondingly) the distance between neighbouring resonant Alfven layers for the 
most interesting region near the MSW resonant point {tmsw ~ O.3i? ). One can 
see that the distance between layers n and n + 1 is changing in order of magnitude 



from z s {n) - z s {n + 1) ~ 200 km (B = 10 G) up to z s (n) - z s (n + 1) ~ 2000 km 
(B = 100 G). 



(42) 



Az. 



s 




(43) 
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5 Density fluctuations 



Considering whole solar interior < z < R & for the same background magnetic field 
values B = 10 G and B = 100 G we have found many resonant layers shown in our 
main Figures 8a, 8b where an unnormalized density fluctuation 5p(z)/p (calculated 
from Eq. flTT|)) is plotted as the function of z, < z < R Q . 

Notice that a height of spikes occurs the random one because of uncertainty of 
numerical calculation near each resonant layer when the numerical code is running 
along whole solar radius with many spikes skipping occasionally some sharp maxima. 

Obviously from comparison of Figs. 8a, 8b that in the case of the small magnetic 
field Bq = 10 G spikes are placed more tightly than in the case of Bq = 100 G 
and really resemble a density matter noise since small distance between neighbouring 
resonant layers shown in Fig. 7a,b. 

Let us estimate the enhancement of the density perturbation Sp/p given by Eq. 

©, 

A = $p/po)z=z s _ I v z (£)/H + u 2 (£) \ z=Zls . s 

(Sp/ Po ) z \v z (0/H + u 2 (0\ z ' 1 ] 

in the vicinity of an Alfven resonant layer, z — z\ s <C H, obeying the inequalities 

(a - 1)< 4 = %e z/H r x « i, 

c 

where the real part of the ^-argument equals to £i = e^ z ~ zls ^ H (here £i > 1 since 
z > z\ s ) and the first inequality means a finite small displacement from the resonant 
point z\ s provided the large parameter (3q = 4 x 10 4 3> 1 for the small magnetic field 
value Bq = 10 G is substituted. 

For simplicity let us also consider the longitudinal wave propagation k \\ Bq, k y = 

0, when the index p is real, or for hydrogen with 7 = 5/3 we find p = 2^J (7 — lj/7 ~ 
1.27. 

One can easily check from Eq. (|1(]) that for conditions chosen above the longitu- 
dinal speed component is of the order 

w-WO+W 1 ) > (45) 

where the singular solution tt2(0 in Eq. (|23|) takes of the asymptotic form «2(£) — 
C 2 J-^) « C 2 (u^i^l)-^/Y(p + 1). 

Thus, the amplification Eq. (HI) of a MAG wave density amplitude Sp/po is given 

by 

i(z - z ls (n))/H) 2 + 4rf 2 (n)! i+^/4 



A~ ^- 1SV ' v ' (46) 

" L 4rf 2 (n) J ' K ' 

or enhancement is really big because of transparency for MAG eigen waves shown in 
Fig. 2, d(n) = Q 2 /Q\ 1, especially for large harmonics n> 1. 
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In Figs. 9a,b we have shown the profile Sp n (z)/po of the selected mode n = 300 
built numerically in the idealized case Bq = 10 G and k y = for which the resonant 
layer position occurs just under the top of the convective zone, z s /R Q = 1 — zu/Rq = 
0.9985 (in our model we neglect turbulent properties of medium in that region). 

One can easily understand the continuous exponential increase of the wave n = 300 
propagating from the center z = R Q up the resonant Alfven layer z\ a ~ 0.01 near 
the photosphere (see Fig. 9b). This follows from the solution Eq. (^) extended on 
whole solar interior and rewritten as U\{z) = e~ z ^ 2H (Aioe lz ^^ H + A2oe~ lz ^^ H ). 

For this largest cavity the wave length of MAG waves along z-axis equals to 

~ z s /300 ~ 2300 km and the enhancement seen from Fig. 9a (for unnormalized 
amplitude) is of a few orders of magnitude ( A ~ 10 2 — 10 3 for the distance X z ~ 
(z — zia)) in good agreement with the analytic estimate Eq. (f46|). 

^From the right side of that resonant layer one can see in Fig. 9a a steep decrease 
of density perturbation described by Eq. (|TTD through the solution Eq. (|2"B] ) vanishing 
exponentially to the surface z = 0. 

6 Discussion 

We have studied MHD (or more correct name: MAG) wave origin of the solar 
matter density fluctuations Sp/p that could be important both for Solar Physics and 
for the MSW solution to SNP. 

We have found MAG-wave eigen modes for cavities within solar interior with 
the spectrum given by the dispersion equations Eq. (|40|) . The lower boundary of 
cavities is the common point - center of the Sun, z = R & , where the longitudinal 
perturbations are absent, v z = b z = 0, and the top (upper) boundary is given by the 
Alfven resonant layer position, z = z s (n), that depends on the node number n, the 
background magnetic field value B and the wave numbers k x , k y . 

A finite amplitude of the density perturbation Sp/p in each Alfven resonant layer 
is determined by the mechanism like the collisionless damping of MAG waves that 
is similar with the Landau damping of Alfven waves = absorption by electrons at 
the (7erenkov resonance u = k x VA = kv e cos9, v e 3> v a ~ 1 cm/s, 9 ~ ir/2. Thus, 
electrons moving almost across Alfvenic layer can absorb these waves. 

Crucial and open questions for future exploration are: 

a) How to enhance MAG eigen modes existing at a negligible fluctuation level 
within cavities deep under bottom of the convective zone and where is the energy 
source for an amplification of MAG waves in solar interior? What is the normalization 
of MAG eigen modes? Full set of ortonormalized functions (like in SHM) is still not 
built for these MAG waves. 

b) Does the dispersion relation that could be obtained via a solution of kinetic 
equations in inhomogeneous magnetised plasma coincide with the dispersion found 
from Eq. (f40|) in the MHD boundary problem? 
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c) Can we apply the simplest analysis given here to another geometry of the 
background magnetic field Bq, how to generalize the analytic approach developed in 
present work in the case of an non- uniform magnetic field Bq(z) and how to take off 
the isothermal approximation T = const! Latter step would allow us to construct 
the sound speed profile c s (z) connected with the inverse problem like in SHM. 

As the necessary step for a generalization and for the checking of the simple 
asymptotic solution found here an independent numerical solution of the same initial 
Heun equation Eq. (O) or even of the whole MHD equation system Eq. (|7|) is needed 
to confirm, first, the presence of spikes shown in our Figs. 8a, 8b. 

An answer on the first question relies on the excitation of a contour through 
the connection (via a non-linear interaction) with the second (neighbouring) contour 
excited by the given generator. In the oscillation theory some mechanisms for energy 
transfer in two- or n-contour systems are well-known. They are based on solutions of 
the telegraph equation in electrodynamics for two (or n-cell) connected contours or of 
the non-linear oscillator equations for two (n-cell) connected pendulums in mechanics. 

The key question whether the energy transfer time would be shorter than the 
collisionless damping time of eigen modes in neighbouring contours has a promising 
solution for high harmonics n 3> 1 for which the temporal damping occurs small in 
our model, d{n) <C 1. 

The non-linear interaction in hydrodynamics is given by the vVv-teim. that was 
omitted in the linear theory. This term can lead to an amplification of the mode 
v\ if the second mode i>2 (excited in a contour displaced in space with respect to 
the first one) was generated by an external source. On this way one can imagine a 
turbulent source in the convective zone (with a spectrum including low frequencies in 
our model ((a few days) -1 ) with the continuous transfer of energy down the bottom 
of the convective zone. 

One can not exclude also another energy source connected with radiative losses 
and due to some spatial inhomogeneities of thermonuclear reactions in central region 
or a diffusion there (for instance due to 3 He diffusion near r ~ O.2i? ). In that case a 
source Q= the rate of all energy density sources and losses arises in the r.h.s. of the 
energy conservation law Eq. (|5|). 

To judge of whether MHD regular waves may be a plausible source of matter noise 
let us note how they can be converted to a noise. Even infinite sum over eigen- values 
(radial numbers n = 0, 1.. and angular moments /) does not't lead to the noise as it 
is argued in |22j for g-modes in SHM. 

The presence of regular waves with a small wave length is only a motivation but not 
realization of the white noise. Really an instability (like buoyance force disappearance 
at the bottom of the convective zone as a source of acoustic turbulence there) can lead 
to the white noise with continuous spectrum due cascade or inverse cascade processes 
with any change of the turbulence scale. This is the non-linear process with the decay 
(or fusion) of regular modes (of number 1 and number 2) obeying the conservation 
laws: uj = uji + u>2, k = ki + k2- In hydrodynamics or in MHD this corresponds to 
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the non-linear term (vV)v omitted in linear approximation in SHM (B$ = 0) and in 
linear MHD. 

Thus, without any assumption about a neutrino magnetic moment we have pro- 
posed a way to modify the MSW solution to SNP The registration of 1-10 days 
variations of solar neutrino flux in underground detectors with large statistics of neu- 
trino events seems to be a signature of MHD wave influence. The search of Day/Night 
(D/N) variations due to Earth effects for neutrino oscillations has been already car- 
ried out in the SuperKamiokande (SK) experiment to check the MSW scenario. It 
seems to be not difficult to recognize a few days periods from the Fouirer analysis of 
the SK data if MAG wave density perturbations would occur large enough (!?). 
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Figure Captions 

Fig. 1 The integration region < £1 = Re £ < £ioo for all three asymptotic 
solutions (downwards- to the center of the Sun). Here = 1/2 is the mathing point 
for the solutions 1*2 (0 and u 3 (£) and £ = 3/2 is the matching point for tti(£) and 

Fig. 2. The ratio of imaginary and real parts of the eigen frequencies d —\ Q2 | 
/Qi in dependence on the node number n 

Fig. 3a and b. The normalized phase velocity v ph /c s = va/c s in dependence on 
the node number n: a for the longitudinal wave propagation k y = and for different 
values of the ambient magnetic field Bq (in Gauss), b for the fixed Bq = 10 G and 
for different angles of the wave propagation k±/k x = 1/ cos a = 1,2,3,5. The low 
(dashed) and high (solid) frequency branches matched at n max are plotted 

Fig. 4. The spectrum / = Ui(n, B , k x )/2rc for the longitudinal wave propagation 
(k y = 0). The ambient magnetic field B = 10 G is chosen. The low frequency 
branches (dashed lines) and the high frequency ones (solid lines) matched at n max = 
2780 are plotted for different node numbers n. For short wave lengths \ x = 2ir/k x the 
characteristic periods of MAG waves T = f' 1 = (1 -=-10 /2Hz)- 1 = 11.57-=- 1.157 days 
correspond to the node numbers n = 300 2000 

Fig. 5a and b. The Alfven layer position z s (n): a for the longitudinal wave 
propagation k y = in dependence on node numbers n for the ambient magnetic field 
values B = 10 100 G, b for the fixed background field B = 10 G in dependence on 
the wave propagation angles, (cos a)" 1 = k±/k x = 1,2,3,5. Dashed curves respond 
to too low frequency branches T = f^ 1 > T beyond the approximation of the Sun 
at rest, T <C T Q = 27 days 

Fig. 6. The width Az s (n) of Alfven resonant layers for the longitudinal wave 
propagation k y = in dependence on the node number n for the ambient magnetic 
field values B Q = 10, 20, 100 G 

Fig. 7a and b. The positions z s (n) for Alfven layers near the MSW resonant 
region zmsw ~ 0.3i? Q for the ambient magnetic field: a Bq = 10 G (the distance 
between neibouring layers equals to z s (n — 1) — z s (n) ~ 200 km), h B = 100 G (the 
distance between neibouring layers equals to z s (n — 1) — z s (n) ps 2000 km) 

Fig. 8a and b. The density profile 8p(z)/p (z) in the case of longitudinal wave 
propagation k y — 0. The spikes (with unnormalized amplitudes) respond to Alfven 
resonant layers with the node numbers: a for the ambient magnetic field Bq = 10 G 
from n ~ 300 near the photosphere z s = R Q up n = 2680 at the MSW region 
z s = O.3i? , b for the ambient (large-scale) field B = 100 G from n ~ tens at the 
photosphere z s = R Q up n = 268 at the MSW resonant point zmsw — 0.3i2 Q . Small 
amplitudes of regular MAG waves that fill each cavity from the center of the Sun 
up the corresponding right boundary z s (n) (up a spike) are not seen between spikes 
while really the wave length within an inner cavity, A 2 = 2n/k z = z s /n, is very short 
since n ^> 1 

Fig. 9a and b. The density perturbation profile 5p n (z)/p for the wave mode 
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n = 300 in the background field B — 10 G: a small-scale structure of the density 
perturbation mode n = 300 near the resonant layer position z s /R Q = 0.9985, b the 
same density perturbation mode within the largest cavity with the upper bound near 
the photosphere, z s /R Q ~ 0.99 (not seen on the large scale < z < R G ), 
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Fig. 1. The integration region < £i = .Re £ < ^ioo for all three asymptotic 
solutions (downwards- to the center of the Sun). Here £i = 1/2 is the matching 
point for the solutions «2(0 and 1*3 (£) and £ = 3/2 is the matching point for 
•ui(0 and w 2 (0- 
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Fig. 2. The ratio of imaginary and real parts of the eigen frequencies d = f^/^i 
in dependence on the node number n. 




Fig. 3a. The normalized phase velocity V p h /c s in dependence on the node 
number n for the longitudinal wave propagation k y = and for different values 
of the ambient magnetic field Bo (in Gauss). The low (dashed) and high (solid) 
frequency branches matched at n max are plotted. 
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Fig. 3b. The normalized phase velocity V p h/c s in dependence on the node 
number n for the fixed Bo = 10 G and for different angles of the wave propaga- 
tion ki_/k x = 1/cosa = 1, 2, 3, 5. 
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Fig. 4. The spectrum / = wi(n, Bo, k x )/2ir for the longitudinal wave prop- 
agation (k y = 0). The ambient magnetic field Bo = 10 G is chosen. The 
low frequency branches (dashed lines) and the high frequency ones (solid lines) 
matched at n max = 2780 are plotted for different node numbers n. For 
short wave lengths = 2ir/k x the characteristic periods of MAG waves 
T = / _1 = (1 -i- 10 fjtHz) -1 = 11.57 -f- 1.157 days correspond to the node 
numbers n = 300 2000. 
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Fig. 5a. The Alfven resonant layer position z s (n) for the longitudinal wave 
propagation k y = in dependence on node numbers n for the ambient magnetic 
field values Bq = 10 -j- 100 G. Dashed curves respond to too low frequency 
branches T = / _1 > Tq beyond the approximation of the Sun at rest, T 
T @ = 27 days. 
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Fig. 5b. The Alfven resonant layer position z s (n) for the fixed background 
field Bo = 10 G in dependence on the wave propagation angles, (cos a) -1 = 
^l/^x — 1? 2, 3, 5. 




Fig. 6. The width Az s (n) of Alfven resonant layers for the longitudinal wave 
propagation k y = in dependence on the node number n for the ambient 
magnetic field values Bq = 10, 20, 100 G. 
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Fig. 7a. The positions z s (n) for Alfven resonant layers near the MSW resonant 
region zmsw & O.SRq for the ambient magnetic field Bo = 10 G. The distance 
between neibouring layers equals to z s (n — 1) — z s (n) ~ 200 km. 
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Fig. 7b. The positions z s (n) for Alfven resonant layers near the MSW resonant 
region zmsw & O.SRq for the ambient magnetic field Bo = 100 G. The distance 
between neibouring layers equals to z s (n — 1) — z s (n) 2000 km. 
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Fig. 8a. The density profile Sp(z)/ po(z) in the case of longitudinal wave prop- 
agation k y = 0. The spikes (with unnormalized amplitudes) respond to Alfven 
resonant layers with the node numbers for the ambient magnetic field Bo = 10 G 
from n ~ 300 near the surface z s = Rq upto n = 2680 at the MSW region 
z s = O.3i?0. Small amplitudes of regular MAG waves that fill each cavity from 
the center of the Sun upto the corresponding right boundary z s (n) (upto a spike) 
are not seen between spikes while really the wave length within an inner cavity, 
X z = 2ir/k z = z s /n, is very short since n>l. 
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Fig. 8b. The density profile Sp(z)/ po(z) in the case of longitudinal wave 
propagation k y = 0. The spikes (with unnormalized amplitudes) respond to 
Alfven resonant layers with the node numbers for the ambient (large-scale) field 
Bo = 100 G from n ~ tens at the surface z s = Rq upto n = 268 at the MSW 
resonant point zmsw = O.3i?0. 
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Fig. 9a. The density perturbation profile Sp„ (z) / po for the wave mode n = 300 
in the background field Bo = 10 G within the largest cavity with the upper 
bound near the surface, z s /Rq ~ 0.99 (not seen on the large scale < z < -R©). 
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Fig. 9b. Small-scale structure of the same density perturbation mode (Fig. 9a) 
near the resonant layer position z s /R® = 0.9985. 



